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Abstract: We use the AdS/CFT correspondence to study the transport properties of mas- 
sive N = 2 hypermultiplet fields in an N = 4 SU(N C ) super- Yang-Mills theory plasma in the 
large N c , large 't Hooft coupling limit, and in the presence of a baryon number chemical po- 
tential and external electric and magnetic fields. In particular, we compute the flavor fields' 
contribution to the stress-energy tensor. We find infrared divergences in the stress-energy 
tensor, arising from the flavor fields' constant rate of energy and momentum loss. We regu- 
late these divergences and extract the energy and momentum loss rates from the divergent 
terms. We also check our result in various limits in which the divergences are absent. The 
supergravity dual is a system of D7-branes, with a particular configuration of worldvolume 
fields, probing an AdS-Schwarzschild background. The supergravity calculation amounts to 
computing the stress-energy tensor of the D7-branes. 
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1. Introduction 

The conductivity tensor aij measures the electrical response of a conducting medium to 
externally applied fields. It is defined by 

(Jj) = Gij Ej 
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where E are externally applied electric fields and (Jj) are the electrical currents induced 
in the medium. Similarly, the thermoelectric conductivity tensor measures the thermal 
response. It is defined as 

(Qi) = OLij Ej 

where (Qi) are heat currents induced in the medium, 

(Qi) = (T t i )- t x(J i ), 

where (T 1 -) are the components of the stress-energy tensor, so that (T^) are momentum 
densities, and fi is the chemical potential, governs the Peltier, Seebeck and Nernst effects. 

In this paper we will use the the anti-de Sitter / Conformal Field Theory correspondence 
(AdS/CFT) [1-3] to study the heat currents that arise in response to external fields in a 
strongly-coupled system. As with other calculations of transport phenomena via AdS/CFT 
(for example the calculation of viscosity in ref. [4]) the over- arching objective is to develop 
solvable toy models that may provide insight into the real-time behavior of strongly-correlated 
systems. 

High-T c superconductors tend to exhibit an unusually large Nernst effect even outside 
the superconducting phase, which may give some hint as to the nature of the still-mysterious 
pairing mechanism at work in these materials. While an AdS/CFT description of high-T c 
superconductivity has not yet been found, AdS/CFT techniques can clarify what sort of 
theories could account for a large Nernst signal. For example, one potential explanation was 
suggested in ref. [5]: a quantum critical point. The physics of systems near their quantum 
critical points is described by a strongly-coupled CFT — a class of theories for which AdS/CFT 
provides many examples. Based on this logic, an AdS/CFT calculation of thermo-electric 
response in a strongly-coupled (2+l)-dimensional CFT was performed in ref. [5]. 

We will focus on a (3+l)-dimensional N = 4 supersymmetric SU(N C ) Yang-Mills (SYM) 
theory plasma at temperature T. We work in the limits N c — > oo with the 't Hooft coupling 
A = g\r M N c fixed, and with A > 1. We will introduce a number Nf of massive M = 2 
supersymmetric hypermultiplets transforming in the fundamental representation of the gauge 
group (i.e. flavor fields). We fix Nf such that Nf <C N c , and work to leading order in Nf/N c . 
With mass-degenerate flavor fields, the theory has a global U(Nf) symmetry, whose U(1)b 
subgroup we identify as baryon number (hence the subscript). We will study the theory with 
a finite U(1)b chemical potential and with external electric and magnetic fields that couple 
to anything carrying U(1)b charge. 

AdS/CFT equates the SYM theory in the limits above with supergravity on the ten- 
dimensional spacetime AdS^ x S 5 , where AdS$ is (4+l)-dimensional anti-de Sitter space and 
S 5 is a five-sphere [1]. The SYM theory in thermal equilibrium is dual to supergravity on an 
AdS-Schwarzschild spacetime, where the SYM theory temperature T is identified with the 
Hawking temperature of the AdS-Schwarzschild black hole [6,7]. The Nf hypermultiplets 
appear in the supergravity description as a number Nf of D7-branes embedded in the AdS- 
Schwarzschild background [8]. With Nf <C N c D7-branes, we may treat the branes as probes 
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and neglect their back-reaction on the geometry. We will explain the supergravity description 
of the hypermultiplet mass, the U (1) b chemical potential, and the electric and magnetic fields 
in the sequel. Though we focus on this system, our analysis easily extends to other probe 
D-brane systems. 

One drawback of the system studied in ref. [5] is that it has translation invariance, which 
implies momentum conservation. The system thus has no way to dissipate momentum, and so 
the DC transport behavior was singular. For example, the DC conductivity at finite density 
is infinite because the charge carriers, in the presence of an external electric field but without 
frictional forces, accelerate forever. In ref. [5], this problem was addressed by introducing 
(by hand) a relaxation time r, representing the presence of defects. A calculation of r in 
an AdS/CFT toy model (a not-quite-consistent truncation of M-theory) was later performed 
in ref. [9]. A calculation of the Nernst effect entirely within the framework of a "complete" 
AdS/CFT system remains to be done. 

The D7-brane probe of ref. [8] provides an elegant solution to the problem of singular 
DC transport in a translationally- invariant finite density system [10,11]. The key fact in the 
SYM theory is that the free energy of the flavor fields and the adjoint fields have different 
scalings with powers of N c . The flavor fields have order NfN c degrees of freedom, whereas 
the adjoint fields have order N%. As an external electric field accelerates the flavor charge 
carriers, their momenta increases at a rate of order N c . In a stationary state (as constructed in 
refs. [10,11]), this external force is balanced by a dissipative force which transfers momentum 
to the adjoint degrees of freedom of the SYM plasma at a rate of order N c . As the energy 
and momentum densities of the adjoint degrees of freedom are of order N%, they can absorb 
the order N c contribution without experiencing any significant back-reaction up to times of 
order N c . That is, the adjoint degrees of freedom effectively act as a heat sink into which 
the flavor fields can dissipate energy and momentum at a constant rate. Only at very late 
times (parametrically large in N c ), will the adjoint degrees of freedom acquire an order-one 
velocity. Before that time, their motion is negligible. As a result, the DC conductivity of the 
flavor fields is finite even at finite density [10, 11]. 

As for the thermo-electric response, we will see below that the issue of dissipation becomes 
more subtle than it was for the conductivity. In particular, no meaningful separation exists 
between the flavor fields' and the adjoint fields' contributions to the heat currents. We 
must face the fact that if external work is done on the system then the heat current is not 
stationary. Nevertheless, we are able to extract many interesting pieces of information from 
the momentum and heat currents. For generic background fields, we calculate the rate of 
change of the heat and momentum currents (which we, irrespective of the sign, will refer to 
as "loss rates") and find they agree with field theory expectations. We also study special 
limits, or carefully crafted observables, that are insensitive to the loss rates. 

The general strategy is to match conserved charges in the bulk of AdS with those in 
the SYM theory. Energy and momentum densities in the bulk map one-to-one to the cor- 
responding densities in the SYM theory. This strategy has been implicit in much recent 
work involving fundamental strings in the bulk which are thought to describe states with 
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large energy or angular momentum, starting with ref. [12]. Here we generalize this idea to 
probe D-branes. We find that in general this strategy sometimes misses important contri- 
butions to the energy of the state which come from the back-reaction of the probe on the 
geometry, and in particular the change in the horizon of the AdS-Schwarzschild black hole 
due to the probe. Based on these observations, we can clarify a few important points in the 
literature. Momentum currents, on the other hand, can reliably be calculated in the probe 
approximation . 

This paper is organized as follows. In section |2] we review the basic strategy and inves- 
tigate the issue of back-reaction with a few illustrative examples. In the remainder of the 
paper we apply the strategy to the probe D7-brane. In section [3|, we review the setup of the 
probe D7-brane in the AdS-Schwarzschild background, describing a SYM theory with finite 
U(1)b chemical potential and background electric and magnetic fields, discuss the polariza- 
tion tensor of the flavor fields, and then compute the stress-energy tensor for the D7-brane. In 
section ||] we use our result for the stress-energy tensor to calculate the energy and momentum 
loss rates, finding perfect agreement with field theory expectations based on hydrodynam- 
ics. Unfortunately, these loss rates also prevent us from obtaining new transport coefficients, 
such as otij. In section || we consider a few special cases in which the loss rates vanish. In 
particular, we show that the zero-temperature limit of our result agrees with field theory ex- 
pectations, giving a nice check of our method. In section ^ we conclude with some discussion 
and suggestions for future research. We collect some technical results in two appendices. 

2. Conserved Charges in AdS/CFT 
2.1 The Hamiltonian Framework 

Witten, in his seminal paper on the basics of the AdS/CFT correspondence [3], provided 
two alternative points of view on how to relate field theory and bulk quantities. The most 
commonly used dictionary is the operator /field mapping, where every operator of the field 
theory is associated to a field in the bulk supergravity theory and the source strength and 
vacuum expectation value of the operator can be read off from the leading and sub-leading 
near-boundary behavior of the bulk field. 

The alternative point of view is the Hamiltonian framework. If the bulk and boundary 
theories are truly equivalent as physical theories, they must have the same Hilbert space, 
and therefore every state in the field theory must map to a state of the bulk theory. As a 
consequence, for any given state, the expectation values of conserved charges such as energy 
and momentum must agree between the bulk theory and the boundary theory. 

This point of view is particularly useful when studying the correspondence in the presence 
of additional sources in the bulk, such as probe D-branes or semi-classical Nambu-Goto strings 
describing single quarks (or quark/anti-quark pairs). In these cases, using the standard 
operator /field dictionary to calculate the total energy requires solving for the back-reaction 
of the probe, so that the stress-energy tensor of the field theory, T^, can then be read off from 
the asymptotic behavior of the bulk metric [13]. Alternatively, the Hamiltonian framework 
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guarantees that we can simply calculate the total energy of the bulk source and equate it 
with the total energy of the boundary solution. 

This strategy has been used successfully for the study of Nambu-Goto strings, in par- 
ticular for spinning strings, following the work of ref. [12]. Here we will apply the same 
philosophy to probe D-branes. We will study spatially homogeneous solutions, so both bulk 
and boundary charges are simply the volume of space times a constant density, and we can 
go one step further and directly equate the densities, rather than just the charges. 

2.2 Back-reaction 

At finite temperature, one subtlety that complicates the analysis is that the black hole horizon 
can carry part of the bulk charge. For example, for an electric test charge moving in the 
presence of a Reissner-Nordstrom black hole, the total boundary charge would be the sum of 
the test charge and the black hole charge. For sources that cross the horizon, as our probe 
D-branes will, the probe makes two contributions to the bulk energy: the energy of the part of 
the probe outside the horizon plus the change in mass of the black hole due to the presence of 
the probe. That is, for the purpose of calculating the energy of the probe, we cannot neglect 
the back-reaction of the probe. 

Consider an action of the general form 

S = N^Sgrav + ToSprobe, (2-1) 

where we have extracted Newton's constant from the gravitational action and converted it 
to the field theory quantity N% and similarly we have extracted an overall normalization 
To from the probe action, where To <C N% (in the case of a string To is the string tension, 
which is order one in the large- N c counting, whereas for a probe D-brane To is of order N c ). 
The back-reaction of the probe will only change the metric by a small correction of order 
e = Tq/N^. Naively, we would therefore expect that the back-reaction is totally negligible to 
leading order in e. However, when calculating the energy density the order e correction to 
the order N% energy density from the gravity part is of the same order, To, as the leading 
contribution from the probe action, so both must be included. 

The only thermodynamic quantity that can be reliably calculated in the probe approxi- 
mation is the free energy, which is simply minus the on-shell action. In this case, the order 
e correction to the order N% on-shell action of the gravity sector is proportional to 5 3 g ™J 
evaluated on the background (uncorrected) metric, which vanishes by the equation of mo- 
tion 1 . The only contribution to the free energy at order Tq is thus the probe's on-shell action 



1 One subtlety in this argument is the potential contribution of boundary terms, which was analyzed in 
detail in ref. [14]. The variation of the Einstein-Hilbert term 16 * G \f§R m the action vanishes in the bulk but 
leads to a finite boundary term of the form 16 * G J dv y/hn^v^ where dV is the boundary of spacetime, n M the 
outward-pointing normal, W the induced metric on the boundary, h its determinant and 

n M u M = n M /i pCT (V p Sg^a - V M Sg pa ). 

The second term above cancels the variation of the Gibbons-Hawking term f g \/hK j '(8ttG). In the standard 
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(evaluated on the background (uncorrected) metric). This fact has been used in most recent 
thermodynamic studies of flavor physics using probe D-branes, for example refs. [15-18], but 
has been less appreciated in the study of string probes. We will illustrate this point in two 
examples below. 

An important conclusion to draw for probe D-brane systems is that for quantities where 
the order N 2 term from the background vanishes (so that the contribution of the plasma dual 
to the bulk degrees of freedom is only of order one) the only contribution at order N c comes 
from the probe D-branes, and we can safely neglect the back-reaction. In particular, this 
allows us to calculate momentum currents in the static black hole background reliably from 
the probe D-brane stress-energy tensor. Energy densities, on the other hand, will receive 
contributions from both the probe and the horizon. 

2.3 Examples 

2.3.1 A Simple Toy Model 

Let us first demonstrate the basic point in a simple exactly-solvable toy model. Consider pure 
Einstein gravity with a negative cosmological constant coupled to a spacetime-filling probe 
D-brane of tension -j^q] for a probe D-brane To /G is a small parameter of order Nf/N c . The 
Euclidean action is 

s -IeW^4* + T-4 (2 - 2) 

In the absence of the brane (To = 0) the equations of motion that follow from this action 
allow planar AdS-Schwarzschild solutions 

ds 2 = h{r) eft* + _ + _ ^ h{r) = __ _A_ (2 3) 

t is the time coordinate and dx 2 is the metric of three-dimensional Euclidean space, while r 
is a radial coordinate that runs from the horizon at r = to the boundary at r = oo. The 
Ricci tensor, Ricci scalar, and temperature associated with this geometry are 

r> __± n p_ d(d + l) T _\_rj L d_ 

L 2 L 2 ' (3~ L 2 Att : ( ' 

so that L is the radius of curvature of the space 2 . 

variational problem determining Einstein's equations one takes the variation Sg^v to vanish at the boundary 
and so the first term which involves its derivative along the boundary vanishes, too. In this case the variation 
of the action vanishes identically and the variational problem is well-defined. In our case Sg^v is the change in 
the metric due to the backreaction of the brane, which in general will not vanish for spacetime-filling branes 
that extend out to the boundary. However, as this surviving boundary term is proportional to a gradient of 
the backreaction along the boundary it will vanish for field theory configurations without spatial gradients, as 
we consider in this work. 

2 This is the only subsection of the paper in which we write the L explicitly. In all other parts of the paper 
we will use units in which L = 1. Also, this subsection and the next are the only places in the paper where we 
use a Euclidean-signature metric. In all other subsections we use a Lorentzian-signature metric. 
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We identify the free energy density of the field theory as minus the on-shell action (divided 
by (3 and the volume J dx), which is 3 



f " "levrGL^+i r " h ~ ~\tc (t) Ld ~ lTd - (2 ' 5) 

The entropy density s, energy density e, and pressure P can be obtained from / via standard 
thermodynamic relations or directly from the gravity side: the entropy is the Bekenstein- 
Hawking entropy associated with the horizon area and the stress-energy tensor was calculated 
in ref. [19]. The results are 

P = -f, e = (d-l)P, sT = dP (2.6) 

as appropriate for a conformal field theory. 

The effects of the tension To can easily be incorporated analytically, as To simply shifts 
the cosmological constant. All thermodynamic quantities only change by replacing the old 
curvature radius L with the new curvature I subject to the relation 

d(d-l) d(d-l) 

T? ~ To = (2 ' 7) 

or, expanding to leading order in t = TqL 2 , 



Correspondingly, the shift in free energy is 

It = ^f. (2.9) 

From Sf we can also calculate Se = — (d — 1)5 f and T5s = —dSf. 

We want to confirm that the effect of To can also be correctly captured directly by a probe 
calculation. As is well-documented in the literature, the probe D-brane's contribution to the 
free energy should simply be given by — fi^ 1 times the on-shell action of the D-brane in the 
original, uncorrected AdS black hole background. Naively, the gravitational on-shell action 
evaluated on the back-reacted solution (that is, the solution incorporating the order t effects 
due to the D-brane) should also contribute at order t, but this contribution is proportional to 
the variation of the action with respect to the metric evaluated on the zeroth-order solution, 



3 Here we have added counterterms to cancel an r — » co divergence of the on-shell action (interpreted as 
a UV divergence in the field theory). The counterterms produce a finite contribution that is — | times the 
contribution from the r — rn boundary of the integral, so that the overall answer is 1/2 times the r = rn 
term. In ref. [6], where background subtraction was used instead of counterterms, the factor of 1/2 arises 
from making sure that both black hole and the thermal AdS background subtraction are really at the same 
temperature. 
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which vanishes because the zeroth-order solution extremizes the gravitational action. Indeed, 
we find that the D-brane's on-shell action gives 



Sfprobe ~ 327rG ° L d-l j ~ 2d f ( 2 - 10 ) 

in perfect agreement with the analytic answer expanded to order t. Calculating the D-brane's 
energy from its stress-energy tensor, which we will denote U^ u , we find 

Se P robe = - f dr^U\ = J drj-g = 5f. (2.11) 

The energy density is equal to minus the on-shell action, instead of the expected (d— 1) times 
the on-shell action. 

Note what is missing: we haven't accounted for the part of the probe D-brane hidden 
behind the horizon. The only effect of that part of the D-brane is to increase the horizon 
area. The total energy in the bulk has to include this change in black hole mass to count the 
total energy in the bulk correctly. Indeed, the first law of black hole thermodynamics, which 
is purely a statement of classical gravity, tells us that the change in horizon area comes with 
a change in black hole mass given by 

Se bh = T5s bh = -dSf (2.12) 

so that indeed 

5e bh + 5e probe = -{d- 1)5 f. (2.13) 

The contribution of the D-brane alone only accounts for a small fraction of the total bulk 
energy, and moreover even contributes with the "wrong" sign. We will encounter this issue 
in various forms throughout this paper. 

2.3.2 Finite- Temperature Correction to Quark Mass 

Our second example, more relevant to finite-temperature calculations performed in the AdS 
dual of strongly-coupled N = 4 SYM, is the mass shift of a heavy quark in the thermal N = 4 
SYM plasma, which was first considered in ref. [20]. As we will explain below, dynamical 
flavors in the SYM theory are represented by probe D-branes in AdS that extend down to 
some radial position r m (in the coordinates of eq. (^1|)). In the gravitational theory, a single 
quark is represented by a string stretching all the way from the horizon to the probe D-brane. 
At zero temperature the energy of such a string is given by m = r m- For this static 
solution, m is identified as the mass of the quark. As we argued above, the reliable quantity 
to evaluate is the free energy, i.e. minus the on-shell action. For the string worldsheet under 
consideration, g rr gtt = 1, so even at finite temperature we find that the on-shell action is the 
length of the string times the circumference of the Euclidean time circle; the only difference 
at finite T is that the string only stretches from the horizon at to r m . The result is that 
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we find for the free energy of the quark 4 

F = -^- i {r m -r h ) = m-]-V\T (2.14) 
lira I 

leading to a temperature-independent entropy of S = Using thermodynamics, we can 

compute the energy of the quark at finite temperature, 

E = F + TS = m. (2.15) 

That is, the energy of the quark is independent of temperature. 

In ref. [20] the energy of the stretched string was calculated directly via the canonical 
momenta of the Nambu-Goto action, with the result E = m — ^y/\T. This is clearly not 
the correct answer for the energy of the string, as it contradicts basic thermodynamics. 
Apparently the back-reaction of the string on the black hole will have to change the mass of 
the black hole by \\f\ to account for the correct energy E of the quark. In ref. [20], \V\ 
was somewhat loosely referred to as the change in mass of the quark. A precise statement of 
what this quantity really represents is the change in free energy of the quark. 

Let us summarize the important lessons of this section. In both of our examples, the 
probe's stress-energy tensor (or, for the string probe, canonical momenta) gave the free energy 
precisely, not the energy. In other words, the probe's stress-energy tensor completely misses 
the entropy contribution. This is to be expected, as the contribution to the entropy takes the 
form of an expanding horizon, which is a result of the back-reaction of the probe. The free 
energy is given by the probe action, the entropy is given by the back-reaction, and the energy 
is given by both the probe action and the back-reaction. 

3. The Stress-Energy Tensor of Flavor Fields 
3.1 The D7-brane Solution 

In this subsection we present the supergravity solution describing massive hypermultiplets 
propagating through an J\f = 4 SYM plasma with finite U(1)b chemical potential and in 
the presence of external electric and magnetic fields. We will also provide a rough sketch of 
the SYM theory phase diagram, to give a sense of where the results of later sections will be 
valid. Most of this subsection is a summary of refs. [10,11]. For readers already familiar with 
refs. [10, 11], we recommend skipping to the next subsection and consulting this subsection 
as needed. 

We use a Fefferman-Graham [21] form for the AdS§ metric 

dz 2 

ds Ads 5 = + 9tt dt 2 + g xx dx 2 (3.1) 

4 The expression below is the answer for a quark of fixed color. For a quark of unspecified color F receives 
an additional trivial contribution of log(iV c ). 
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where z is the AdS$ radial coordinate. The boundary is at z = 0. When we need an explicit 
form of the metric, we will use 



*U = " ^\/Jf£ *? + > + zA/4) (3 ' 2) 



The black hole horizon is at z = zh with zjj 1 = T. Here we are using units in which 
the radius of AdS§ is equal to one. In these units, we convert from string theory to SYM 
quantities using ol~ 2 = A. 

We will use an S 5 metric of the form 

ds 2 s5 = dd 2 + sin 2 9 ds 2 sl + cos 2 9 ds 2 s3 (3.3) 

where 9 is an angle between zero and it/ 2 and ds 2 sl and ds| 3 are metrics for a unit-radius 
circle and 3-sphere, respectively. 

We next introduce Nf probe D7-branes. The relevant part of their action will be the 
Born-Infeld term, 

Sm = -N f T D7 J cPCy/-det (dab + (2na')F ab ). (3.4) 

Here is the D7-brane tension, C, a are the worldvolume coordinates, g a b is the induced 
worldvolume metric, and F^b is the U(l) worldvolume field strength. The D7-branes will be 
extended along all of the AdS§ directions, as well as the S 3 directions inside the S 5 . 

Our ansatz for the worldvolume fields will include the worldvolume scalar 9{z). The 
D7-brane induced metric is then identical to the background metric, except for the radial 
component, which is g zz = \ + 9'(z) 2 , where prime denotes differentiation with respect to 
z. The supergravity field 9{z) is dual to the SYM theory operator which, roughly speaking, 
is the mass operator of the flavor fields. We will denote it as O m . The precise form of O m 
appears for example in ref. [22]; just thinking of O m as the mass operator will be sufficient for 
our purposes. We can identify the mass and expectation value (O m ) from 9(z)'s asymptotic 
form, 

9{z) = c lZ + c 3 z 3 + 0{z 5 ). (3.5) 

The mass m of the flavor fields is then m = 7^7, and the expectation value is given by 
(O m ) oc-2c 3 + icf [17,23]. 

We discuss 0(z)'s equation of motion and boundary conditions below, around eq. fl3.19|) . 
We will soon discuss other D7-brane worldvolume fields, in particular worldvolume electric 
and magnetic fields; 9(z)'s equation of motion will of course depend on their values. We can 
in principle solve for 9{z) numerically, but in what follows we will not do so. Instead we 
will write our results in terms of 9{z) and make use of special cases where we know exact 
solutions. For massless flavors {c\ = 0), 6(z) = is always an exact solution for any values 
of the external fields, and corresponds to massless quarks with (O m ) = 0. Our D7-brane 
will have a worldvolume magnetic field, however, so we must remember that, as shown in 
refs. [24-27], in the presence of a background magnetic field 9 = is typically only meta- 
stable, and the true ground state has a non-trivial 9{z) even when c\ = 0. In the field theory 
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we interpret this as spontaneous chiral symmetry breaking induced by the magnetic field. 
Still, we will use the meta-stable 9 = solution as a simple example to confirm that our bulk 
expressions agree with field theory expectations. 

In the AdS-Schwarzschild background, two topologically distinct types of D7-brane em- 
beddings exist, the so-called "Minkowski" and "black hole" embeddings. In a Minkowski 
embedding, 9{z) describes a D7-brane for which the S 3 inside the S 5 "slips" as the D7-brane 
extends from the boundary into the bulk of AdS§, eventually collapsing to zero volume, i.e. 
9{z) begins at zero at the boundary and reaches ^ at some z > zh, so that the D7-brane 
smoothly terminates outside the horizon. In a black hole embedding, 6{z) never reaches ^ 
outside the horizon, so the D7-brane intersects the horizon. In the absence of worldvolume 
electric and magnetic fields, a first-order phase transition from Minkowski to black hole em- 
beddings occurs as m/T is lowered through a critical value [15-18,28,29]. The transition is 
often called a "meson melting" transition: Minkowski embeddings give rise to a gapped, dis- 
crete meson spectrum in the SYM theory, while black hole embeddings give rise to a gapless, 
continuous meson spectrum [30]. 

The U(Nf) gauge invariance of the coincident D7-branes is dual to the U(Nf) symmetry 
of the mass-degenerate flavor fields in the SYM theory. We identify the £7(1) subgroup as 
baryon number, U(1)b- The D7-brane worldvolume Abelian gauge field is dual to the 
SYM U(1)b current J M , so to introduce a finite U(1)b density in the SYM theory, we must 
introduce the worldvolume gauge field A t (z). We may identify the U(1)b chemical potential 
[i and the density (J 4 ) from the asymptotic form 5 of At(z): 



where the constant c gives the density: (J*) oc c [22]. 

As shown in ref. [22], with nonzero At(z), only black hole embeddings are allowed. The 
physical reason is simple: with nonzero At(z) the D7-brane has an electric field pointing in the 
radial direction, F z t- These field lines are sourced by a density of strings. The force that such 
a density of strings exerts on the D7-brane overcomes the tension of the D7-brane, pulling the 
D7-brane into the horizon. The strings are then "hidden" behind the horizon, where the field 
lines may safely terminate. Numerical [22] and analytic [31] analysis of D7-brane embeddings 
has confirmed that with nontrivial At(z) the only physical D7-brane embeddings are black 
hole embeddings. When we work with nonzero (J 4 ) in the SYM theory, our D7-brane will 
always intersect the horizon. 

The first-order "meson melting" phase transition discovered at zero density has been 
shown to persist to finite density, producing a line of first-order transitions that ends in a 
critical point [22]. A region of the phase diagram near the line of transitions is also known to 
be thermodynamically unstable [22]. Our results will not apply in the unstable region. 




(3.6) 



5 A good question is why /i is physical: it is constant, so normally we could gauge it away. We discuss the 
boundary condition on At (z), which forbids such gauge transformations, as well as the boundary conditions 
on the other gauge fields, in appendix A.]. 



To introduce perpendicular electric and magnetic fields, and the resulting currents (J x ) 
and (J y ), we also include in our ansatz the gauge field components [10,11] 

A x (t, z) = -Et + h x (z), A y (x,z) = Bx + h y (z). (3.7) 

In each case, the leading term is a non-normalizable mode that introduces an external field 
into the SYM theory. Choosing a gauge in which A z = 0, we can write the nonzero elements 
of F ab as: 

F tx = -E, F xy = B (3.8) 
F zt = A' t , F ZX = A' X , F zy = A' y . (3.9) 

To date, the SYM theory in the presence of external fields has been analyzed (via 
AdS/CFT) with only E [27,32] or B [24-27] alone, not both simultaneously, and only with 
zero U(1)b density. As mentioned above, at zero temperature the magnetic field triggers 
spontaneous chiral symmetry breaking and also produces a Zeeman-like effect on the me- 
son spectrum. The magnetic field also increases the meson melting transition temperature. 
Heuristically, the magnetic field "holds mesons together," and exhibits a kind of "magnetic 
catalysis" of chiral symmetry breaking. Indeed, in supergravity language, for sufficiently large 
magnetic field only Minkowski embeddings exist. On the other hand, an electric field lowers 
the meson melting temperature, which is easy to understand intuitively: the electric field will 
pull quarks and anti-quarks (or squarks and anti-squarks) in opposite directions, so that we 
do not need to heat up the system as much to "melt" mesons. 

Our ansatz only involves gauge field components in the four-dimensional (z, t, x, y) sub- 
space, so the D7-brane action takes the form of the four-dimensional Born-Infeld action, times 
some "extra" factors coming from the extra dimensions, 

S D7 = -N j dz cos 3 9 g l J 2 ^j-g - (2ira') 2 ^gF 2 - (27ra') 4 ^ (F A F) 2 . (3.10) 

We have divided both sides of eq. ( |3.10| ) by the volume of M 3,1 , so that technically Sm 
is an action density. We have defined g = g zz gu g xx as the determinant of the induced 
metric in the (z, t, x, y) subspace. We have also introduced the constant M, which, using 

T m _ a '~ i 9s 1 _ XNc ;„ 
lD 7 ~ ~l2^p 2^F' 

M ee NfTmVs* = jAiN f N c (3.11) 

where Vgs = 2ir 2 is the volume of a unit-radius S 3 . Writing F 2 = F^F^, where Greek 
indices run over (z,t,x,y), we have explicitly 

2 9 F 2 = glx A 't + 9tt9xxA x + gttgxxA' 2 + g zz g X xE 2 + g zz g tt B 2 (3.12a) 
-(FA F) 2 = B 2 A' 2 + E 2 A 12 - 2EBA[A' y . (3.12b) 
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Starting now, L will denote the Lagrangian density, albeit with an unconventional sign choice, 
introduced for future convenience: Sd7 = — / dzL. 

Clearly the action only depends on the z derivatives of At, A x and A y , so classically the 
system has three constants of motion. As shown in refs. [10,11], we can identify these as the 
components of the U(1)b current density in the SYM theory, 

(f) = §. (3.13) 

Our system thus has a nonzero U(1)b density (J 1 ) as well as U(1)b currents (J x ) and (J y )- 
Given these constants of motion, we can solve algebraically for the derivatives of the gauge 
field (the field strength components): 



A , iz) = _ V^\(jy-Ba K{z) = x m ff>| (3.14) 
9xx J&T^ V \9tt\ ^X^ 

V m\ Vtx - a 2 

where 

£ = \9tt\9L + (27Ta') 2 (\g tt \B 2 - g xx E 2 ), a = (2rra') 2 (\g a \(Jt)B + g xx (J y )E) (3.16) 

X = \gtt\g 3 xx M\^a') A cos 6 9 + (2na') 2 (\g tt \(J t ) 2 - g xx ((J x ) 2 + (J y ) 2 )). (3.17) 

In terms of these quantities the on-shell action, which can be identified as minus the free 
energy density in the grand canonical ensemble, is given by 

S D 7 = -AA 2 (2vrQ') 2 J dz cos 6 6 g 2 xx ^fg zz \g tt \ (3 ' 18) 

The equation of motion for the worldvolume scalar 9(z) can be derived in two ways, either 
from the Legendre-transform of the action, which we denote as <Sd7, 

= " (2^)2 / dz 9ll 2 \9tt\~ 1/2 9xxV£x- a 2 , 

or by varying eq. ( |3. 10j ) and then plugging in the solutions in eq. (ETUI) and ( pUD . 

To complete the D7-brane solution, we must specify boundary conditions on 9{z). For 
Minkowski embeddings, we have 9(z m ) = ^ for some < z m < zh and 9'{z m ) = 00 to avoid 
a conical singularity [17]. For black hole embeddings, the boundary conditions are 9'{zh) = 
for the embedding to be static, while 0(zh) G [0, is a free parameter. Notice that each 
case has a free parameter (the value of z m or 9(zh)) which maps in a one-to-one fashion 
to the asymptotic parameter c\ oc m. For example, with nonzero At(z) (hence only black 
hole embeddings), the limits are 9{zh) = maps to c\ = (when B is nonzero, this is the 
meta-stable massless solution) and 9{zh) — ► f maps to c\ — > 00. 
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3.2 The Polarization Tensor 

In the next subsection we will use the D7-brane solution above to compute the stress-energy 
tensor of the flavor fields. Some contributions to the stress-energy tensor come just from the 
electric polarization and the magnetization of the medium, however, as we will now explain. 

Even in an equilibrium system, background electric and magnetic fields produce non- 
vanishing momentum currents due to polarization effects. As reviewed for example in ref. [5] 
(and references therein), even in equilibrium we expect a contribution to {T^ u ) of the form 

(T^) pol = M^Fl. (3.20) 

where M^ v is the polarization tensor, 

= (3.21) 

with f2 the free energy density (and where we take the derivative with other variables held 
fixed). The components of M^ a with one t index and one spatial index are electric po- 
larizations while components with two spatial indices are magnetizations. The full energy- 
momentum tensor (Tt) then divides into two pieces: 

to = to fimd + m PO i (3.22) 

where, for example, (T^) /hu<i corresponds to the genuine momentum current due to the flow 
in the medium. As reviewed in ref. [5], both (Tt) and (Tt) fluid obey the same conservation 
equation, 

d»{T llv ) = F vp {J< } ), (3.23) 

but only (Tt) fluid represents observable quantities that can couple to external probes of the 
system, and hence is the appropriate object to use when studying transport, for example, 
when computing transport coefficients. In particular, we should use (T^) to identify the 
heat current densities (that we discussed in the introduction), 

{Qi) = {T\) fluid ~li{Ji)- (3.24) 

In gauge-gravity duality, we identify = —Sd7, where here Sd7 is the D7-brane action 
evaluated on a particular solution for the worldvolume fields. For our choice of background 
electric and magnetic fields, the x and y components of the polarization, M u (i = x,y), and 
the magnetization M xy will be non-vanishing: 

M* = M xy = (3.25) 

dEi dB 

Notice that both the electric and magnetic fields must be nonzero for the polarizations to 
contribute to the momentum densities: 



{T\) = {T l x )f luid - M ty B, (T* y ) = (T^) fi u id + M tx B. 



(3.26) 
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Notice also that M ty will be nonzero even though we have not introduced a background 
electric field in the y direction. Suppose we did introduce an electric field in the y direction, 
E y . We then easily find that, due to the F A F term in the action, eq. ( [3.10D , taking the 
derivative with respect to E y and then setting E y = produces a nonzero result. Indeed, a 
little algebra shows that 

dS D1 _ 1_ 
dE v ~ E 



dz (A' x (J y ) ~ A' y (J x )) 



(3.27) 



which we will use to simplify expressions in the sequel. A useful fact to remember is that 
M ty = when (J*) = 0: when (J*) = we can use the result for the conductivity in appendix 
[A] to show that (J y ) = 0, and from the explicit solution for A' y in eq. ( 3.15| ) we can also show 
that (J*) = implies A' y = 0. 

The calculation of M tx and M xy from S&7 is more complicated. Consider for example 
M tx . In field theory terms, we need to compute jg\ T B - We start with eq. ( 3.10|) , evaluated 
on a particular solution. The on-shell action Sd7 = — ^ will then have explicit E dependence, 
as well as implicit dependence through the solutions for 9(z) and the worldvolume gauge 
fields. We thus employ the chain rule 6 , 



dS D7 
dE 



dz 



dL 89 dL 89' 8L ^ 



=t,x,y 



8A\ 8L 
8E 8 A'; 



(3.28) 



Notice in paricular that in the ^ term the derivative only acts on the explicit E dependence 
in L (i.e. on the explicit factors of E appearing in eqs. (|3.12a|) and (|3.12b|) ). We then use 
the fact that partial derivatives commute to write = -§^-§ei an d integrate by parts to 

find 



dS D7 
dE 



dz 



dL 
8E 



+ 



89 8L 
~8E~89' 



■-U 



8L_ 
~89 ~ 

E 

i=t,x,y 



8 8L 
~8z89' 

8 A, 8L 



89 
~8E 



E 



8Ai 8 8L 
~d~E~dzdN i 



8E 8 A', 



(3.29) 



Obviously, of the terms under the integral, the term in parentheses and the terms in the sum 
over i vanish due to the equations of motion. That leaves the ^ term under the integral, 
and the boundary terms. We will not compute ^ in what follows, nor will we discuss the 
boundary terms. Our main point is that the only contribution to the polarization from the 
bulk of AdS§ comes from jM. A similar statement applies for the magnetization M xy , 



i.e. 



the only bulk term comes from jjg. In sections |5.3| and |hj we will take limits in which the 
boundary terms vanish (or are negligibly small), in which case the bulk terms become the 
only contributions to the polarization tensor. 

In the next subsection we will see factors of -§jjj-, and ^ appearing in the stress- 
energy tensor. Most of these arise from the expected contribution to (T^ u } from (T^ u ) po i. 

6 We are using arguments similar to those in refs. [18,26,33]. 
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3.3 The Stress-Energy Tensor 

Our goal now is to compute the contribution that the flavor fields make to the stress-energy 
tensor of the SYM theory, using the above holographic setup. We have chosen to work in the 
Hamiltonian framework. If pi denotes the momentum associated with the flavor fields in the 
SYM theory, with i = x,y our momenta of interest, then in the Hamiltonian framework we 
identify the conserved charges 



Pi 



J dtdx (T*) = J dtdxdzd 3 a yf^g^ U\. (3.30) 



The a are coordinates on the S 3 wrapped by the D7-branes and gD7 is the determinant of 
the induced metric on the D7-branes. We have introduced the notation U t i as the D7-branes' 
momentum density, reserving the notation (T^) for the expectation value of the flavor fields' 
momentum density in the SYM theory. Note that in order to form a covariant quantity, we 
must include a factor of yJ—QDi i n the integral. In principle, a similar factor should also 
appear in the four-dimensional integral, but as the SYM theory lives in flat space, this factor 
is unity. Furthermore, if the energy-momentum tensors are independent of the four spacetime 
coordinates, then the integrals over dtdx will only produce a factor of the spacetime volume, 
so that we can equate the momentum densities directly: 



(T\) = j dzd 3 a yf^gl^ U\. (3.31) 



Our task is thus to compute the stress-energy tensor of the D7-branes. Let us first define 
notation. Let 

Q a b = J dzd 3 a ^g~^U a b . (3.32) 

When the indices a and b are in SYM theory directions, we can identify {T a b ) = Q a b . The 
indices a and b can also be in the z or S 3 directions, however, in which case the SYM theory 
interpretation requires more effort. In what follows we will be able to translate some, but not 
all, such components into SYM theory quantities. 

We can compute Q a b in two different ways. We can of course directly compute the 
variation of the D7-brane action, Sd7, with respect to the background metric. Alternatively, 
because the momenta are the generators of translation symmetries, we can derive the tensor 
components via a Noether procedure. The two methods must agree up to boundary terms 7 . 
Indeed, we have used both methods and have found perfect agreement (not just agreement 
up to boundary terms). As the calculation by variation of the action is somewhat lengthy, 
we include it in appendix |B|. The result of the Noether procedure is 

e- l --/*(L* + 2F**L-«^), (3.33) 

where we have performed the trivial integration over the S 3 , since our ansatz for the world- 
volume fields is independent of these directions. 



7 See for example sections 7.3 and 7.4 of ref. [34]. 
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Note that the bulk theory does not actually have translation invariance in the z direction. 
Calculating T^ z using the Noether procedure may thus seem suspicious. As mentioned above, 
however, even for these components we have verified explicitly that the Noether calculation 
agrees with the calculation via variation of the background metric. 



We expect the last term in eq. (|3.33[) to contribute to T z z , given our ansatz 6{z). We 



find, however, that the last term in eq. ( 3.33| ) also contributes to the T l z components with 



i = t,x, y. To see why, suppose we allow 9 to depend on t, x, y. We then find that, due 
to the gF 2 term in eq. ( |3.10 ), taking the derivatives Jjfg, j^f-Q, and -^-g and then setting 



dtO = d x 6 = d y 6 = produces a nonzero result. This is very similar to what we saw in the 
last subsection, where the y polarization was nonzero even though our ansatz has E y = 0. 
We write explicit expressions for the derivatives -^-g in appendix ^. 

We will now present all the components of the stress-energy tensor. 

Many of the components are simple. For example, in the S 3 directions, and in the third 
Euclidean field theory direction, the only components are on the diagonal, and all are simply 
— J dzL = Sd7- 

The nontrivial components are in the (z,t,x,y) subspace. Notice that, with one index 
up and one down, the energy-momentum tensor will not be symmetric, so we computed all 
sixteen components separately. We will also identify current components, (J*), (J x ) and (J y ), 
whenever possible to simply the expressions. 

We will now write all of the components of Q a b in the (z, t, x, y) subspace. For notational 
simplicity, in what follows, we will not write the J dz, which appears for every component. 
Primes will denote Jg. 

The components with upper index t are 

Qt t = -L- F xt ^- F zt -§^ = -L + #f§ + ( J*)^ 
9', = -F„jfc - F yx ^ = {,P)A' X - | ((J*)A> y - {JV)A' X ) 

c*t _rp SL _ rp SL p dL i / Tt\ a l 

^ y ~ r %y&F tx z V8F tz ~ D dE^\ J i J± y 

r\t 77i SL rp SL I Qi SL _ a I dL ai dL . gi SL 

V z — r ^z SFtx r yzsF ty 5d t e — ^xdE ^ydEy Sd t 8 

The components with upper index x are 

Q x t = ~ F *inb = ( JX ) A 't 

0X T 771 SL 77 SL 77 SL T I TP dL _i_ p dL I / TX \ A I 

x - -L - ^txsFT t - t V^jFT v ~ b **wz ~ ~ L + tj m + B m + W ) A x 

C\x 77 SL / jx\ A 1 

w y ~ rz vSF xz - W l A y 

0x _ p SL _ rp SL _i_ qi SL At dL , ai dL , qi SL 

z — r tz x tp . r vz xtp 1 " xa a — »tp ~r jt-,, an T 



tz SF xt r y z SF xy ~r u sd x e ~ A t 8e A y dB ~r u sa x e 
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The components with upper index y are 



= -Fzxwt ~ Ftx M~t ~ ( JX ) A 'y 



@\= -L-F xy ^-F zy ^ = -L + B% + (jy)A' y 

r\V _ _ T7i SL _ T7i SL i gi SL _ a' dL ai dL_ , pi SL 

yj z — r tz & F yt rxz 5F yx i " t7 5d y 8 — A tdE y dB Sdyd 

The components with upper index z are 

® Z *= -Ftojfa-Fyzjfc = {J*)E + {JV)B 

@z z = -L- Y,i=t,x, y F izm- + e 'w = ~ L + T,i=t,x, y ( Jl ) A 'i + e 'w 

All quantities on the right-hand sides are evaluated on-shell. Those components with both 
indices in field theory directions (t, x, y) we can identify with the energy- momentum densities 
of the flavor fields, as explained above. 

We would like to convert the components of Q a b to field theory quantities. In most cases, 
whether we can do so depends on whether we can perform the z integration. Sometimes this 
is easy. For example, we know that J dz L = —Spj = £1, and J dz A' t (z) = —fi, where \i is 
the U(1)b chemical potential. We thus have, for example, @ x t = (T x t ) = —fj,(J x ). In some 
cases we can translate to SYM theory quantities without doing the z integrals. For instance, 
terms with §^ or ^ multiplying E or B we interpret as contributions from the polarization 
tensor, as explained in the last subsection (and as will be verified explicitly, in certain limits, 
in section ||). In the next subsection we will identify the & z i components (i = t, x, y) with the 
rates of energy and momentum loss of the flavor fields. Notice also that the Q z z component 
is, up to the O'tet term, identical to the Legendre transform in eq. ( 3.18| ). 



On the other hand, we have not found a field theory interpretation for the components 
@ z , Q x z , and Q y z , for which the z integration is non-trivial. For many components (® t x , 
Q t y , Q x x , @ x y , etc.), converting to SYM theory quantities requires integrating A' x and A' , 
and the field theory meaning is not immediately clear. We discuss the z integration of A' x 



and A' y in appendix A.l 



Finally, notice that @ t x = when (J*) = 0, partly because 8*^. includes the expected 
polarization term BM ty , and as explained in the last subsection, M ty = when (J*) = 0. 
That Q t x vanishes when ( J 1 ) = is easy to understand physically. As explained in appendix 
[A], our system has two types of charge carriers, the charge carriers we introduced explicitly 
in the density (J*), but also charge carriers produced by pair production in the external 
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electric field E. When (J*) = 0, we find a nonzero charge current {J x ), coming from pair 
production. However, we expect the momentum current in the x direction to vanish because 
charges produced in pairs will have zero net momentum. 

4. Energy and Momentum Loss Rates 

As soon as we turn on any external fields in the plasma, the work done by the external 
forces will change the energy and momentum of our system at a constant rate. Without a 
mechanism for dissipation, a stationary solution is impossible, as the momentum (and energy) 
in the system will never stop increasing. The plasma is translation-invariant, so momentum 
can never really dissipate. However, in the probe limit Nj <C N c , the flavor degrees of freedom 
are very dilute relative to the much more abundant J\f = 4 SYM degrees of freedom. The 
current, momentum and energy densities of the flavor fields can be constant, but they will 
transfer energy and momentum density to the N = 4 plasma, which will thus gain energy 
and momentum at a constant rate. The conservation law for the full stress tensor is (see for 
example ref. [5]) 

d^(T^)=F up (J"). (4.1) 

In particular, for the spatially homogeneous solutions we considered, with only F xt = E and 
F xy = B nonzero, we have 

d t {T\) = -E(r) (4.2) 
d t (T t x ) = E(J t ) + B(jy) 

d t (T%) = -B(r). 

From the point of view of the D7-branes' stress-energy tensor, these loss rates are reflected 
in the appearance of an IR divergence, very similar to what happened for the dragging string 
solution of refs. [20,35]. To be more specific, Q t t , Q t x and Q t y have divergences from the 
z = zh endpoint of integration (in addition to any expected UV divergences from the z = 
endpoint, which can be cancelled with counterterms) . 

If we regulate these IR divergences by only including those parts of spacetime that had 
time to communicate with the boundary within a time t, we find that the divergences are 
linear in t and the coefficients can be interpreted as loss rates. Let us demonstrate this 
quantitatively. We simply need to expand the expressions from the previous section and 
study the behavior close to the horizon. Let z = z^ — e. We find that in all three cases the 
integrand has terms that diverge as K The loss rate should be proportional to the coefficient 
of that pole, but we must fix the exact proportionality constant. To do so, we can compare 
the i poles to the time required for a ray of light to propagate from z = z^ — e out to the 
boundary, which is given by 
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The integrand in tn g ht itself diverges as ^ producing the \ pole. We can thus extract finite 
rates by identifying the divergent parts of Q t i as tug^t times the rate. The rates so obtained 



are in perfect agreement with eq. (4.1). 



The right-hand sides in eq. ( |4.2| ) are closely related to Q z i with i = t,x, y. Specifically 
the Q z i differ only by a factor of the z integration, Jq H dz = zh = y , from the loss rates on 
the right-hand sides in eq. ()4.2|). 

In ref. [36] , the response of the SYM theory plasma to a moving source, such as a heavy 
quark, was computed using the usual AdS/CFT framework, that is, by introducing a probe 
source in AdS$, such as a long string, computing the back-reaction on the metric to linear 
order in the perturbation caused by the probe, and then extracting the SYM theory stress- 
energy tensor from the asymptotic form of the back-reacted metric. One of the main results of 
ref. [36] was that the non-conservation of the stress-energy tensor (the right-hand sides of eq. 
( [4.2| )) was determined by the u zi" components of the probe's stress-energy tensor. We have 
reached the same conclusion using the Hamiltonian AdS / CFT framework, which confirms the 
equivalence of the two approaches. 

Given the above loss rates, we can construct two IR-safe quantities, 

h = E(T\) - B{T\) 

h = (J*)(T* X ) + (JV)(T* V ) + <J*><T* 4 ). (4.4) 

Both I\ and I2 are free of IR divergences, so the quantities they represent on the field theory 
side are constant as long as the currents are time-independent. We can write these two IR-safe 
objects as the t components of the currents 

1$ = e x ^(T»)F pa (4.5) 
1% = (Tt){r) (4.6) 

where X3 refers to the third Euclidean spatial direction in the field theory. We can easily show 
that these two currents are conserved. The divergence of If is proportional to F A F and 
so vanishes in topologically trivial background fields, such as our orthogonal E and B. (An 
interesting generalization of our results, and those of ref. [10,11], would be to the case with 
non-trivial E-B.) Using d^T^) = F v P(J p ), we find that the divergence of 1% is F pu {J^){J v ) , 
which vanishes identically. 

As conservation of 1^ seems to be a peculiarity of the background we chose, we will focus 
on the interpretation of 1% , which is conserved whenever the currents are stationary. Indeed, 
I2 has a simple interpretation: it represents the mass-energy current density 4-vector of the 
fluid as measured by inertial observers with a 4-velocity proportional to (J v ), (see e.g. section 



3 The limit of zero temperature but finite density, with E > B, is more subtle. In this case, the system is 



not really stationary, as we review in more detail in appendix A. 2. Indeed, in that case the loss rate seems 
to acquire more terms. Of course what really happens at zero temperature is that the real physical loss rate 
is no longer time-independent, as the charge carriers continuously accelerate. The loss rates from 0^ on the 
other hand still agree with the hydrodynamic expectation. 
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4.2 of ref. [37]). At least in the limit where the current is only carried by quasi-particles, we 
can understand this result in simple terms. In the rest frame of the charge carriers, the force 
due to the external fields does not lead to any work, as W = J F ■ ds = (the charge carriers 
are at rest and the fluid moves). Such an effect has already been seen in comparing the 
results for the dragging string of refs. [20,35] to the calculation of ref. [38], both describing 
a quark pulled by an electric field and moving at constant speed through the plasma. In the 
rest frame of the fluid, a quark that has been moving for a long time builds up an energy 
proportional to the distance traveled [20,35]. In the rest frame of the quark, the same quark 
has a time-independent, finite energy. 



5. The Limit of Zero Temperature 

To extract transport coefficients from currents we need a stationary system. For generic 
external fields, our system is not stationary, as visible from the loss rates discussed above. In 
certain limits, however, the loss rates may vanish or be higher order in the external fields. In 
this section, we consider such limits, and check that our results agree with the expectations 
of the SYM theory. 

A special case in which all the loss rates vanish is when all the currents vanish, so that 



we are analyzing equilbrium physics. As we discussed in section |3.2| , polarization effects due 
to the external fields may still produce non-vanishing momentum currents even when the 
currents (J x ) and (J y ) are zero. At zero temperature, the two cases where we can set all 
currents to zero consistently are E = and B > E. From the supergravity perspective, these 
are the only cases where (J x ) = (J y ) = and the D7-brane action remains real for all z. 

A third scenario in which the loss rates can be neglected, even for E > B, is when E ~ e 
and B ~ e for some very small e. As e carries dimension what this scaling really means is that 
we can neglect the loss rates on time scales smaller than In this case, the loss rates 

at zero temperature are of order e 5//2 , whereas the leading contribution to the momentum 
density is of order e 2 , hence to leading order the loss rates are negligible. This is due to the 
somewhat peculiar scaling (J x ) ~ E 3 / 2 of the current at zero temperature (see appendix |A.2| ). 
At finite temperature, where the current has typical Ohmic form (linear in the electric field), 
we have not been able to identify an interesting scenario where the loss rate can be neglected. 
All the examples we will present thus involve the specialization to zero temperature. 

5.1 E = 0, B and (J*) nonzero 

As a first example where the loss rates vanish, we simply set E = 0, so that ( J x ) = {J y ) = 0, 
and the system is in equilibrium. In our D7-brane solution, only At(z) and B are nonzero. 
From our result for Q l t = {T l t ), we find, when E = 0, 



Q\ = - J dz (L-(J*)4) (5.1) 

which, using J dz L = —Sd7 = ^, immediately implies 

n = -(T t t )- f i(J t ). (5.2) 
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At zero temperature we do not expect a horizon contribution to (T* t ), hence (T\) should 



accurately represent the energy density. Eq. (5.2) is indeed the correct thermodynamic 
relation between the energy density -(!r t ) and the grand canonical free energy density. 

Notice that our general result for Q t t remains the same at finite T. At finite temperature 
eq. (^^) is not the right relation between free energy and energy density, however. We 
are missing a Ts contribution, where s in the entropy density, for which we would have to 
calculate the change in horizon area due to the back-reaction. This is in perfect analogy 



to what we found in the toy model of section |2.3.1f the stress-energy tensor of the probe 
completely misses the entropy contribution. 



5.2E = 0,B = 0, (J 1 ) and (J x ) nonzero 



At zero temperature, if all external fields vanish then the system has no dissipation. We may 
of course introduce (Jt), but also constant (J x ) and (J y ). At zero temperature the system 
is boost-invariant, so we can always boost to a rest frame in which only (Jt) is non-zero. 
Without the electric field, nothing distinguishes x from y, so for simplicity we only present 
the result for (T^ x ), since the expression for (T* y ) is identical. 

At zero temperature we have g xx = \gu\, while g zz = -4j + (0 1 ) 2 is generically nontrivial. 
Fortunately, as we show in appendix [B| (see eq. ( B.16| )), (T*,) does not depend on g zz : 



(T t x ) = -(J x ) [ dzA' t ^- = —(J x ) I dzA> = (J x }». (5.3) 
J \9tt\ J 

With this we can determine the heat current, 

{Q x ) = {T t x )-p{J x ) = 0. (5.4) 

This result is actually dictated by Lorentz invariance, as can be easily seen via field 
theory arguments. In the rest frame of the charges, with density ( J*) resi = p, energy density 
-{T\) = e rest , and pressure (T x x ) = (T y y ) = P rest , we know that (J x ) rest = (2*.) = 0. 
Boosting with velocity v and boost parameter 7 _1 = \/T— v 2 , we find 

(J,) = (J x ) = viP, (J 1 ) = IP, (5.5) 

and hence the momentum density in the boosted frame is 

{7* x ) = (e rest + P rest ) V1 2 = (J x )p, (5.6) 

where in the last step we used the equilibrium relation e rest + P rest = p rest p and we needed to 
absorb one power of 7 into p = fi restr y. In a relativistic theory p is best viewed as the constant 
expectation value for the t component of a background gauge field A^. As such it transforms 
nontrivially under Lorentz boosts. That our supergravity result correctly reproduces the SYM 
theory expectation is encouraging. 
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5.3 B > E, (J*) = 



If the density of charge carriers is zero and the magnetic field is larger than the electric field, 
then (J x ) = (J y ) = and the loss rates vanish (see the discussion in appendix |A.2|) . We do 
not expect to find any transport coefficients at zero temperature and density, but we should 
be able to understand the stress-energy tensor in terms of polarizations. In particular, we 
expect (ry = 9^ = M tx B. (Recall that with (J*) = we have Q\ = -M*»B = 0, as 



explained at the end of section 3.3 



To verify that (T y ) = M B, we will focus on the meta-stable 9 = solution representing 
massless flavor fields. Prom our result for the stress-energy tensor with (J*) = we have 
@\ = Bf dz§^. When (J x ) = (JV) = and 6 = 0, the boundary terms in eq. |T2f) vanish, 
so we can identify M tx = J dz ^ and hence ® t y = M tx B. As expected, in this limit only 
the polarization term contributes to the momentum density. 

5.4 E > B but both E and B small, (J*) = 

Finally, if we assume E ^> B, but that both E and B are of the order of some small parameter 
e, and assume again that (J*) = 0, then at zero temperature all loss rates will be of order e 5 / 2 . 
If we work to order e 2 , all quantities will be free of IR divergences and we should again find 
only polarization contributions to the stress-energy tensor. We will also set 6{z) = again, 
for simplicity. 

With (J*) = we have (J y ) = and again Q x = 0. We would like to confirm that 
Q f y = M tx B to order e 2 . Now in eq. ( |3.29| ) the A' x boundary term is non- vanishing, so 



(J* 



dA 



dE 



Zll 



(5.7) 



ft will be proportional to E (the A acts on the E 2 term in eq. ( 3.12a| )). We also have 
IJsHo^ = We Jo** d z A' x , where the leading E dependence of A' x is contained in (J x ) (see 



eq. ( |3.14p ). As explained in appendix A. 2, at zero temperature (J x ) ~ E 3 / 2 so in our limit 



( J x ) ~ e 3 / 2 . The boundary term is thus of order ( J x ) d ^f E ' ~ e 2 , so B ~ e times the boundary 
term is order e 3 , and thus we find 

{T\) =e\ = B j dz^ = M tx B + O (e 3 ) . (5.8) 
Up to order e 2 , only the polarization term contributes to the momentum density. 



6. Conclusion 



We began with the intention of extracting thermoelectric transport coefficients from the heat 
currents of J\f = 4 SYM theory with flavor fields. However, we have found that the probe 
brane approximation does not adequately separate the physics of the flavor charge carriers 
from the physics of the M = 4 SYM plasma. As a result, the loss rates of the charge carriers 
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provide infinite contributions to the heat currents: for this supergravity background, and for 
our particular ansatz for D7-brane worldvolume fields, stationary solutions that would allow 
us to define the time- independent thermoelectric coefficients do not exist. 

Using the natural time-scale of the problem tu g ht, we were able to extract the loss rates 
and compare them to field theory expectations, to which they conformed perfectly Also, we 
were able to calculate the heat currents in specific zero-temperature cases and confirm that 
they had the expected structure in terms of currents and polarizations. 

Another important lesson that we learned along the way is that the probe brane stress- 
energy tensor alone does not give the full order Nf/N c contribution to the energy of the 
system. In order to calculate the energy to this order, the back-reaction of the probes on the 
horizon must be included. 

A possible next step would be to include flavor fields in the fluid dynamical formulation 
of ref. [39]. By looking at the back-reaction in the form of a slowly- varying horizon, the SYM 
theory thermoelectric transport coefficients could be extracted. 
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A. The Conductivity 

In this appendix we will review the results of refs. [10, 11] for the conductivity tensor of the 
flavor fields, as computed using probe D7-branes in AdS/CFT. We then briefly discuss the 
boundary conditions on the D7-brane worldvolume gauge fields, and some subtleties about 
the zero-temperature limit of the result for the conductivity. 

In the SYM theory, we expect (J x ) and (J y ) to be fixed once we choose E, B and (in 
the canonical ensemble) (J*). The main observation of refs. [10,11] was that only one choice 
of (J x ) and (J v ) allows for the gauge fields and the on-shell action to be real for all values 
of z between the horizon and the boundary. More specifically, the three functions £, x an d 
a defined above in eq. d3ltl) all have a zero between the horizon and the boundary, and 



the only way for the action and the gauge fields to remain real is if all three functions pass 
through zero simultaneously. We thus get three equations "for free." The location of the 
zero, which we denote z*, is given by £ = 0. In terms of the dimensionless quantities 

I E \ B 

e = -{2-kol')Ez 2 h = — =— , b = -(2kc!)Bz\ = — =— (A.l) 
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the explicit form of is 



i±- = e 2 - b 2 + yV - 6 2 ) 2 + 2(e 2 + b 2 ) + 1 (A.2) 

4 



-J ((e 2 - 6 2 ) + yV - & 2 ) 2 + 2(e 2 + b 2 ) + l) -1. 

At that point, g 2 ^- 2 *) = 7T T ^"(e, 6) where 

JP( e , fe) = I (l + e 2 _ 6 2 + ^( e 2 _ 52)2 + 2 ( e 2 + 52) + . (A . 3) 
The remaining two equations, a = x = 0, then fix the values of the currents, 

(J x ) = 2 ^ff^^, 2R2 V(gL + (27ra0 2 g 2 )AA 2 (27^a0 4 g^cos6g(^) + (27^a0 2 (J f ) 2 (A.4a) 



(2ro , ) 2 (J i )5 
+ (27ra') 2 S 2 

with all functions of z evaluated at z*. Converting to field theory quantities, we have the 
conductivity, 



'jV*N?T2 ^3/2 _ ^ i p 2^ 



°"~ = V 16vr 2 fr 2 +^ C ° S6 ° {Z * ] + WTW' (A ' 5a) 



= (A ' 5,)) 



where p is defined similarly to e and b, 

, (J*) 



f^T 2 ' 



(A.6) 



Notice that cr^ clearly has two terms adding in quadrature, and that even if we set 
(J*) = (so p = 0), we find a nonzero a xx . Even without the density of charge carriers (J*), 
if we impose an electric field we will see a current ( J x ), from the cos 0(z*) term in a xx . In other 
words, the system has some other charge carriers besides those in (J*). What is the source of 
these other charge carriers? The claim of ref. [10] was that they come from pair production in 
the external electric field. The primary piece of evidence is the behavior of 6{z^). As explained 
at the end of section [O], 9(z#) = maps to m = in the SYM theory, and 6>(z*) — > | maps 
to m — > 00. In the former limit we expect pair production to be maximal, and indeed the 
current coming from pair production is a maximum: cos 0(2*) = 1. In the other limit we 
expect pair production to be suppressed, and indeed we see cos#(z*) — > 0. Notice also that 
the pair-produced charges do not contribute to the Hall current: a xy is proportional to (J*). 

In appendix [A.2j we discuss the zero-temperature limit of the conductivity. 
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A.l Boundary Conditions on Gauge Fields 

The D7-brane worldvolume gauge fields have the following asymptotic behavior: 

a ^ = »-\nB*? z2+ ° {z% (A - 7a) 

A x (z) = -Et + c x + I J^L- Z ^ + 0(z% (A.7b) 

A y (z) =Bx + c y + \ M ^ al)2 * 2 + 0(z 4 ), (A.7c) 

where c x and c y are constants. 

In the usual AdS/CFT recipe, we fix the leading coefficients and then "integrate into 
the bulk," where some boundary condition then fix the solutions completely, and hence fix 
the values of the sub-leading coefficients. This process maps onto the field theory process in 
which, once we choose the parameters in the Lagrangian, the dynamics of the theory fixes 
expectation values of operators. From this point of view, our calculation of the conductivity 
looks strange. We solved for the field strengths (see eqs. ( 3.14| ) and ( |3.15| )) and then imposed 



reality of the D7-brane action to fix the values of (J x ) and (J y ). Implicitly, we must be 
imposing some boundary condition on the gauge fields somewhere in the bulk of AdS§. In 
this section we briefly clarify the boundary conditions on the gauge fields. 

We first consider At. The leading, non-normalizable constant we identify as the U{\)b 
chemical potential, [x. In this case, the background geometry forces a boundary condition 
on At(z): the Killing vector corresponding to time translations becomes degenerate at the 
horizon, hence for the gauge field to remain well-defined as a one-form we must take At(zn) = 
[22]. Notice that we are then also restricted to the subset of gauge transformations that 
leaves this boundary condition invariant, which is why fi has physical meaning (we cannot 
simply gauge it away). 

A x and A y are more subtle. For these, the leading, non-normalizable terms include — Et 
and Bx, which give rise to the gauge-invariant field strengths F tx = E and F xy = B, but the 
leading terms also include the constants c x and c y . We will choose a gauge in which these are 
zero. What then is the IR boundary condition that fixes (J x ) and (J y )? The field strengths 
A' t (z) and A' y {z) are of course gauge- invariant, as are their integrals, 

dzA' x (z)=h x (z H ), / dz A' y (z) = h y (z H ). (A.8) 
o Jo 

Here the background geometry does not force any particular values for h x {zn) and h y (zn) 
upon us. Instead, we are implicitly choosing these constants to produce precisely the values 
of (J x ) and (J y ) that keep the action real for all z. The values of h x (zn) and h y (zn) 
corresponding to these values of (J x ) and (J y ) could be computed explicitly by performing 
the above integrals, using the solutions in eqs. (|3.14 ) and ( |3.15 ) and the values of (J x ) and 



(J y ) in eq. ( A. 4 ). Crucially, however, these integrals diverge in our setup. 
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Many components of the stress-energy tensor that we computed in section || can only be 
translated to field theory quantities using J dzA' x and J dzA' y . As we discussed in section 
|], the divergences in these integrals are related to loss rates in the SYM theory, and indeed 
we expect some components of the stress-energy tensor to be sensitive to the loss rates. The 
conductivity does not depend on the values of f dzA' x and J dzA' y , however: it is insensitive 
to the loss rates. 

A. 2 The Conductivity at Zero Temperature 

Here we study the zero-temperature limit of the conductivity, eq. flA.5| ), in detail. 
The zero-temperature limit of the quantity T defined in eq. ( |A.3 ), 



T(e, b) = - (l + e 2 - b 2 + vV - b 2 ) 2 + 2(e 2 + b 2 ) + lj (A.9) 

is subtle. Recall that an inverse power of T 2 is hidden in the definitions of e and b, so these 
quantities grow large as T becomes small. Correspondingly, the term under the square root 
is dominated by y (e 2 — b 2 ) 2 = \e 2 — b 2 \. Due to the absolute value, we must take care to 
distinguish between the B > E and the B < E cases. 

For B > E, the leading order 1/T 4 piece in T vanishes, and T goes to a constant: 
T = B 2 /(B 2 — E 2 ). In this limit, a xx vanishes and a xy = (Jt)/B as required by Lorentz 
invariance, since for B > E we can boost back to a system with E = 0, hence no spatial 
currents and i?-field only. 

On the other hand, for E > B we find T = e 2 — b 2 . With this both a xx and o xy have 
a finite limit. Specializing for example to the zero- mass, zero-density case of 9 = p = 0, 
we find {J x ) ~ E^l 2 . The system has no linear (Ohmic) current, having instead a current 
that scales as a fractional power of the electric field. Indeed, this power of E is dictated by 
the scale invariance of the zero-density, zero-mass theory. The strong electric field polarizes 
the medium so that even in the absence of charge carriers the medium can conduct at zero 
temperature. 

One note of caution: these finite zero-temperature conductivities can be misleading. For 
example, in the limit of large mass and large density, where a xx is completely dominated by 



the second term under the square root in eq. (|A.5|) , at zero T we find a xx = (Jt)/E, or in 
other words (J x ) = {Jt)- In this case the dynamics is easy to understand: the finite density 
of charge carriers accelerates forever, and hence comes closer and closer to the speed of light. 
Both (J x ) ~ i>7 and (J t ) ~ 7 diverge individually as the boost factor 7 increases. Their ratio, 
however, approaches a constant. In this case the finite conductivity does not represent a truly 
stationary system. The case with zero {Jt), however, seems to be truly stationary. 

B. Stress-Energy Tensor from Variation of the Action 

Our goal in this appendix is to compute the stress-energy tensor of the D7-branes directly, 
by varying the D7-branes' action with respect to the background metric. We do so for two 
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reasons: first, to check our calculation via the Noether-derived form, eq. (|3.33| ), and second, 
to provide explicit formulas that may be useful for future calculations of the back-reaction of 
the D7-branes. 

Let us first fix our notation. In the rest of the paper, g ab denoted the induced metric of 
the D7-brane. In this appendix only, g ab will denote the metric of the background spacetime 



in which the D-brane is embedded. For the particular ansatz in section 3.1, in which the only 
nontrivial worldvolume scalar was 6{z), the distinction was irrelevant except for g zz . The 
induced metric had g zz = + but was otherwise identical to the background metric. The 
distinction between the two is crucial for computing the stress-energy tensor. Physically, we 
want to compute the D7-brane's response to a variation of the background metric, not to a 
variation of its induced metric. Also, in this appendix only, we will denote multiplication by 
2ira' with a tilde. For example, F ab = (27ra')F ab . 

We begin with an expression for the D7-branes' energy-momentum density, 

IT- + **»), (B.l) 

where we have taken a symmetrized derivative with respect to the external metric g ab to 
guarantee a symmetric tensor. Lowering one of the indices in eq. (|B.1|), and plugging into 



eq. ( 3.32 ), we find 

9» c = / dzSa g bc . (B.2) 

J V °9ab og ba J 



™ab = 9cd -^--E^b~ + m = det ( m ab), (B.3) 



Let us introduce some more notation. Let 

dX c dX d 

where the first term in m ab is the induced metric of the D-brane. The X a represent the 
worldvolume scalars and the Q a represent the worldvolume coordinates. We then have simply 

Sdj = -N f T D7 J d\ V^- (B.4) 

(B.5) 



We can then use the chain rule 

8m 8m 8m c d 
8g ab 8m cd 8g ab 

and the identities 

mim'^dc (B.6) 



8m cd dX a dX b 8m 



8g a b d( c dC d ' 8m cd 
to compute the variation with respect to the external metric: 

8S D7 _ N f T D7 1 8m _ N f T D7 m{m- l ) dc dX a dX b 
8g a b 2 J^m~ 8g ab 2 ^m~ d( c d( d ' 



(B.7) 



Our expression for the variation of the D7-brane action, eq. ( B.2 ), becomes 



3 N f T D7 m{m- l ) de f dX a dX b dX b dX a \ 
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which is the main result of this appendix, and does not depend on our particular ansatz for 
the worldvolume fields. Indeed, formally, upon changing the dimension and the tension, this 
formula is valid for any D-branes described by the Born-Infeld action in a background for 
which the dilaton and axion are trivial. A nontrivial NS two-form would simply add a term 
to our definition of m a b- 

To compute @ a b for our particular worldvolume fields, we need to write the matrix m a b 
explicitly and compute the elements of its inverse. To write m a b explicitly, we will need an 
explicit metric for the S 3 , 

0000000 

ds S 3 = da 1 + sin a\ da 2 + sin ai sin a 2 da 3 . (B.9) 

For our ansatz, m a b is 



^ 9tt 


—E 





-4 
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9xx 


B 
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9xx 


-A' 
















A' t 


A' x 
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{g zz + e' 2 ) 
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cos 2 9 sin 2 cti 







V o 

















cos 2 ( 


? sin 2 a\ sin 2 


«2 / 



(B.10) 



which results in 

\/—m = ^g xx cos 3 9 sin 2 a\ sma.2\/d(z), (B-ll) 
where we have defined the shorthand notation 

d(z) = - 5 -(2W) 2 i 5 F 2 -(27ra') 4 i(FAF) 2 

\9tt\gL(9zz + e' 2 ) + (9zz + e ,2 )(B 2 \ gtt \ - E 2 g xx ) (B.12) 

-9 2 XX (A>) 2 + \g tt \g xx ([A' x ) 2 + (A' y ) 2 ) - [A' y E - A! t B) 

where in the first line g is the determinant of the induced D7-brane metric in the (z, t, x, y) 
subspace (i.e. minus the first term in the second line). 

Let us explicitly compute one example, O*^, for which we have 

6* x = I dzd 3 a Nf ^ D7 L- {m{m- l ) xt + m{m- l ) tx ) g xx , (B.13) 
J 2 v-m 

hence we compute 

mim" 1 )^ = g xx cos 6 9{z) sin 4 a x sin 2 a 2 (-g xx A' t A' x - Eg xx (g zz + 9' 2 ) - E(A' y ) 2 + BA' t A'^j 
m{m^) tx = g xx cos 6 9(z) sin 4 a x sin 2 a 2 (-g xx A! t A! x + Eg xx (g zz + 9' 2 ) + E(A'y) 2 - BA' t A' y ^ . 
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With these expressions, and performing the integration over the S s , we find 

Q\ = - f dzM cos 3 9 g 5 J x 2 4^ 

J \MX> 

which, using the definition of (J x ) in eq. ( 3.13] ), can also be written as 

Q\ = -(J x ) J dz A' t 
The same procedure then gives, for example, 



\9tt\ 



(B.15) 



(B.16) 



dz AAcos 3 6g% 2 



{g zz + 8' 2 )BE + g xx A' t A' y 



@\ = - j dz M cos 3 Ogli 2 



9tt\(9zz + 0' 2 )( 9 L + B 2 ) + \g tt \g xx ( (A' x f + (A'f 



Vd(; 



With straightforward algebra, we have checked that all the components of Q a b computed in 
this fashion agree exactly with those derived via the Noether procedure. Checking the Q l z 
components with i = t,x,y requires the following derivatives: 



5L 
5d t 9 


= -AAcos 3 ( 


igli 2 


8L 
5d x 9 


= -AAcos 3 ( 




5L 
5d y 9 


= +M cos 3 ( 


hlL 2 



9xx E A' x , 



v 7 ^ 



(B.17) 
(B.18) 
(B.19) 
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